Developing a robust approach to implementing non-Abelian anyons and topological 
quantum computing in a modified Kitaev honeycomb lattice model 
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Quantum computation provides a unique opportunity to explore new regimes of physical systems 
through the creation of non-trivial quantum states far outside of the classical limit. However, such 
computation is remarkably sensitive to noise and undergoes rapid dephasing in most cases. One 
potential solution to these prosaic concerns is to encode and process the information using topo- 
logical manipulations of so-called anyons, particles in two dimensions with non-Abelian statistics. 
Unfortunately, practical implementation of such a topological system remains far from complete, 
both in terms of physical methods but also in terms of connecting the underlying topological field 
theory with a specific physical model, including the imperfections expected in any realistic device. 
Here we develop a complete picture of such topological quantum computation using a variation of 
the Kitaev honeycomb Hamiltonian as the basis for our approach. We show the robustness of this 
system against noise, confirm the non-Abelian statistics of the quasi-particles to be Ising anyons, 
and develop new techniques for turning topological information into measurable spin quantities. 
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Topological quantum computation shows promise for 
providing an approach to implementing logical quantum 
information operations that are robust against local noise 
sources due to a combination of topology, temperature, 
and a gapped Hamiltonian with appropriate below-gap 
statistics [1-4]. However, open questions remain for 
any physical implementation of such a topological sys- 
tem, as robustness is only guaranteed in the thermody- 
namic limit, and in many cases methods for implementing 
key operations such as preparation and measurement of 
quasi-particles remains at best uncertain 0, . Fur- 
thermore, there does not currently exist an approach for 
connecting potential real-world errors introduced by de- 
fects in Hamiltonian implementation and in operations to 
errors within the topological sector, though preliminary 
efforts to examine errors have commenced 0, |T(| ■ 

In this work we develop a specific physical picture for 
computation and noise using a variant of the celebrated 
Kitaev honeycomb model (llj . Our work follows in the 
spirit of recent attempts to explore explicit architectures 
for topological quantum computation [12j, but here we 
focus instead upon finding appropriate methods for map- 
ping known solutions from circuit-based fault-tolerant 
quantum computing to the case of finite-size, finite-time 
topological computing. In one crucial respect, our work 
departs from prior analyses of the Kitaev honeycomb 
which have focused directly upon the extended version of 
the model, with explicit 3 spin interactions [H; E3 13 > or 
upon small size (16 spin) versions [15| . Instead, we choose 
to introduce a gap in the so-called gapless phase using 
two-spin interactions with ancillary, gap-protecte d sp in 
"gadgets" [16], rather than the original suggestion [Tl| of 
a homogeneous external magnetic field. With this ad- 
dition, the gap introduced is the first perturbative cor- 
rection regardless of the number of vortices (topological 
excitations) existing in the system, in contrast to the case 
of a magnetic field-induced gap. We remark that our ap- 



proach could be implemented experimentally in manner 
similar to Ref. 

We proceed by first developing our approach for im- 
plementing the underlying, extended Kitaev model com- 
prising two-spin interactions only and include all effects 
in the relevant perturbation theory up to the order of 
the induced gap. We then detail how adiabatic control 
of the parameters of the Hamiltonian enables the creation 
and braiding of anyons [17] sufficient for observing Ising 
anyon statistics using numerical approximations based 
upon the exact solution provided first by Kitaev and in 
the spirit of Ref. [18|. We describe an explicit method 
for taking advantage of non-topological sectors, such as 
when vortices are adjacent to each other, to convert topo- 
logical information into spin observables. Finally, we an- 
alyze how control errors can lead to a quantum random 
walk, and show this is suppressed exponentially as vor- 
tex density decreases. Thus we provide a complete pic- 
ture, including errors, for finite-size, finite-time topolog- 
ical quantum computing. 

We start with the Kitaev honeycomb model [ll| , which 
comprises a honeycomb lattice of spin- 1/2 particles with 
two-spin interactions along the links of the honeycomb 
(Fig. [1} . For simplicity of labeling, we deform the hon- 
eycomb into an isometric brickwall structure, such that 
each spin is labeled on a cartesian grid with two spins of 
the unit cell adjacent on the x axis. Thus N x N unit 
cells map to N lines of 2N spins with z — z couplings be- 
tween the lines and alternating x — x and y — y couplings 
along the lines. The corresponding Hamiltonian is 

M N 

m—1 n—1 
+ T rr v rr y 

+ ^ <T 2m-l,n cr 2m,n+l) 



(1) 



One key advantage of this model is the ability to exactly 
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Figure 1: Kitaev honeycomb lattice shown as a tilted brick- 
wall, with each brick corresponding to a single honeycomb. 
Each black and white circle corresponds to a spin, with a 
black-white pair forming the unit cell. Interactions between 
spins are determined by the color of the link (red, x — x; yel- 
low, y - y; green, z - z). 

diagonalize the Hamiltonian by appropriate use of con- 
stants of motion (vortex operators) and use of Majorana 
fermions. 

The Kitaev model is known to have two different topo- 
logical phases, characterized by the relative values of the 
coupling constants. If J a + Jp > J 7 does not hold for all 
permutations of x, y, z, then the model supports Abelian 
anyon excitations and has gapped fermionic modes. Oth- 
erwise it has been suggested that the Kitaev model 
should support non-Abelian vortices but have gapless 
fermionic modes [Tl| . In this paper we will focus on this 
non-Abelian (gapless) phase. 

The anyonic excitations associated with the topology 
can be labeled by defining the vortex operators: W m>n = 

<7 2m-l,n (7 2m,n+l ,T 2rji+l,ii+l (7 2m+2,ti+l' 7 2m+l,ii (7 2ro,ni 

which has eigenvalues ±1 and commute with the Hamil- 
tonian. If the expectation value of W m n is "1; then 
there is a vortex at the plaquette (m, n). For simplicity, 
one can choose J x = J y = J z = J for the gapless phase. 

In the gapless phase, the vortices do not have well- 
defined statistics as they are strongly coupled to the 
fermionic modes near the singularity of the spectrum. 
To open a gap in the gapless phase, and thus enable a 
topologically protected non-Abelian anyonic system, one 
approach is to break the time-reversal symmetry of the 
model [11] by the introduction of a homogeneous mag- 
netic field V h = J2 m ,n M<,„ + + < „).^ We will 
choose a different approach, but first we motivate our 
choice. 

Specifically, there are several often overlooked compli- 
cations to the introduction of the magnetic field. Kitaev 
considered the case where the iV 2 vortex operators had 
eigenvalues of +1, i.e., the vortex-free sector. In this 
case, the introduction of the magnetic field has no effect 
until third order in perturbation theory, leading to the 
introduction of three-spin interactions with a strength 
K ~ h 3 /J 2 , which create a gap of similar size. In con- 
trast, when there are vortices in the system, as would 
occur during any braiding operation, additional terms 
enter at even first order in perturbation theory. This 
can be seen by examination of the situation with two 
vortices separated at a lattice distance much larger than 




Figure 2: Example "gadget" for implementing the effective 
three-spin interaction. 6 ancillary spins are introduced, with 
strong z ~ z (green) and weak x — x (red) interactions. The 
ancillary system is in turn weakly coupled to the underlying 
honeycomb sites with with x — x, x — z, and x — y interac- 
tions (blue, going clockwise). After removal of the high-energy 
subspace of the gadget, the effect is to produce three-spin in- 
teractions between the honeycomb spins. 



the correlation length £ ~ J 3 /h 3 . Within this space, 
of fixed vortex number, the vortices are allowed to hop. 
To be explicit, suppose there is a vortex on the left 
side of the link cfm-i „cr| m n+1 and far away from the 
other vortices (with the state denoted by \tpi)), we know 
that (<Al|Wro-l,n|V>l) = -1 while (V'l|Wm,n|V'l) = +1. 
The magnetic field introduces terms, such as h z a\ m _ x n , 

h*hy x y , hyhyhz x « y / 

J u 2m-l,n°2m,n clIlu JZ u 2m-l,n u 2m,n u 2m+l,n \ LC 

suiting from first, second and third order perturbation re- 
spectively), are of higher or of the same order as the gap 
terms, no matter how small the magnetic field is. E.g., 
the term ft^cfro-i n induces hoping by moving the vor- 
tex to the other side of the link cfm-i n erf m n+1 . Specif- 
ically, since (VM CT 2m-l,n^m-bn CT 2m-l, J^l) = + 1 while 

has a vortex at the plaquette (to, n) instead of (to — l,n). 
This hopping occurs at first order in perturbation theory, 
while the gap is third order. This is particularly challeng- 
ing to work with, even if the gap remains formally cor- 
rect, as in order to manipulate the vortices and to avoid 
unwanted braiding operations, one needs to know and 
fix the positions of the vortices. The vortex dynamics 
induced by the homogeneous magnetic field thus appear 
to go gainst the efforts for manipulating and observing 
non-Abelian anyons experimentally. 



Here we suggest a different approach to generate the 
same three-body terms in the generalized Kitaev model 
using two-body interactions only. Our choice will lead to 
the necessary three-spin interactions to introduce a gap 
in the so-called gapless regime as the first perturbative 
correction that is non-zero, regardless of vortex sector. 
Specifically, we wish to integrate out the auxiliary spins 
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and recover the generalized Kitaev model: 

M N 

H(t) = '5~^X'Jm.n(t) <J 2m.n <7 2m+l,n 
m— 1 n—1 

+ ^ri,ri(^) cr 2m-l,n f7 2m,n 

+ ^m,n(*) (7 2m-l,ii' T 2m,ii+l 

+-^m,nW (T 2m+l,n+l cr 2m,n+l' 7 2m-l,ti 

+-^'m,n(0 <T 2m,n+l f7 2m+2,n+l f7 2m+l,ri+l 

+-^m,n(0 <T 2m+2,n+l cr 2m+l,n+l lT 2m+l,ri 

+-^m,n(0 <T 2m,n cr 2m+l,ri lT 2m+2,ri+l 

+-^m,n(') (7 2ra+l,n (7 L-l,n (7 2m ) n 

+-^ra,n(') (7 2m-l,ii cr 2m,n Cr 2m+l,ii+l) (2) 

We use the idea of perturbative gadgets [2(| to cre- 
ate each 3 spin interaction in Eq. This leads to an 
overhead of 6 additional spins for every 3-spin interaction 
on the honeycomb. Furthermore, we extend the gad- 
get approach by adding an additional protection against 
local noise, thus maintaining the gap- and topological- 
protection. Specifically, rather than using a single set 
of ancillary spins per interaction term, we use two sets 
of auxiliary spins to generate each three-body interaction 
(more sets resulting in better protection). The two sets of 
auxiliary spins are coupled together, as shown in Fig. [3J 

A 

U f Z-Z ~Z i ~Z~Z I ~z~z 

#aux,ifei = - 2"(°j CT fc + a k a l + a l a j 

+ a*a*+a* k a? + at<t*) (3) 

Each set of three spins have ferromagnetic, Ising-like cou- 
plings (z-z). The two sets of spins are coupled to each 
other perturbatively through x-x bonds. The ground 
state manifold of i? a ux,jfci (for A = 0) can be written 
in terms of four 'cat' states, 

i|ttt±iU)®lttt±ttt> (4) 

The x — x bonds breaks the degeneracy of the cat states. 
The states \ |ftt + til) ® Ittt + ttt} and \ |ttt - til) ® 
Ittt - ttt) are still degenerate ground states of the ef- 
fective Hamiltonian, but separated by an energy ~ 3-jp- 
from the other cat states. The two ground states can- 
not be continuously changed into each other under local 
noise without raising the energy, and the transition be- 
tween the two states is suppressed. Furthermore, if such 
noise does cause a change between the two degenerate 
ground states, the effect on the extended Kitaev hamil- 
tonian occurs as local noise, which should be protected 
against if the system exhibits topological order [3|. 

Now we can use gadgets in one of the two degenerate 
ground states to generate three-body terms. We couple 
the original spins to the auxiliary spins by 

Bcoupww = — A(5f of + ofal + of of), (5) 



Suppose each two sets of the auxiliary spins have already 
been initialized to the state 5 Ittt + III) ® Ittt + ttt}- 
Using third order perturbation theory, we obtain 

HeS = -J £ ofofl-l^ E °M°f- (6) 

a-links {j,k,l}£TI 

Compared with Eq. @, we have K = To open a 

large spectral gap comparable to J in the gapless phase 
while preserving the stability of the cat states of auxiliary 
spins, we require J ~ -4? an d A 3 3> A 3 (rather than A 3> 
A). The interaction between vortices can be neglected as 
long as they are separated from each other by a distance 
^> I j^| ~ 1. Finally, control of the three-spin interaction 
is enabled by changing a single interaction between the 
gadget spins and the honeycomb spins. 

In order to study the topological properties of the Ki- 
taev honeycomb model, we generalize it by allowing for 
the coupling terms n , n externally controlled, and 
variable, parameters. The generalized Kitaev model can 
be simplified by spin-fermion transformation (erf = ibfci, 
where bf and Cj are Majorana fermions) which re- 
sults in the Hamiltonian 

H{ f) = lJ2 A jk(t)cjC k , (7) 

where A^ k (t) = —2Jpj* (t)uj k + 2_K£|* (t)v,jiU[ k if j, k are 
linked directly or through a third spin I, otherwise 0. We 
can choose the initial coupling parameters Jp/^O) = Jo, 

and K^(0) = K . The link operator u jk = ///,' /-? 
has eigenvalues ±1. The vortex operator is W p — I!^*:- 

Suppose the system is initialized to the vortex-free 
ground state. We can adiabatically change some of the 
coupling parameters Jp ?£(t) and Kp'*(t). This does not 
change the actual vortex configuration, as vortex opera- 
tors W do not change their expectation value. However, 
it does change the fermion spectrum, and the best de- 
scription is given by considering the same model with- 
out vortices but with a change of sign of appropriate J's 
and K's as suggested in earlier work [18|. These equiv- 
alent descriptions we label pseudo- vortices for this work, 
where the key distinction is that measurement of the vor- 
tex operator does not reveal pseudo- vortices. We remark 
here that the equivalence between changing J and u was 
pointed out in Kitaev's paper to show that the ground 
state energy does not depend on the signs of J [ljj . 

To create a pair of pseudo-vortices at plaquettes 
(mi,n) and (m2,n) in the horizontal direction, we can 
continuously change the sign of all of the n , K^, 6 n 
(mi + 1 < m < m 2 ), and K^ 4 n (mi < m < m 2 — 1) 
between the two pseudo- vortices. We show in Fig. [3] the 
spectrum of lowest energy levels versus the value of n 
for an 18 x 18 lattice (two spins per plaquette, choosing 
mi = 1, m,2 = 10, and n — 10), with periodic boundary 
condition, Jq = 1 and Kq = 0.2. There is a (almost) zero 
energy fermion mode when we have changed the sign of 
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Figure 3: Spectrum of lowest energy levels vs the value of 
Jm,n for creating a pair of vortices on the top right of the 
spins at (1, 10) and (19, 10) in an 18 x 18 lattice (two spins 
per plaquette), with periodic boundary condition, Jo = 1 and 
K = 0.2. 

the relevant parameters between the two vortices. When 
the size of the system and the distance between vortices 
goes to infinity, the first fermion mode should be of ex- 
actly zero energy-in fact, this first fermion mode corre- 
sponds to the topological degeneracy of the two-vortex 
sector. It is also interesting to observe how a nonzero 
fermion mode becomes a zero mode. In fact, all the 
higher modes come down to replace their lower neigh- 
bor modes. 

To move a pseudo-vortex, we can change the coupling 
parameters of the relevant links between the initial po- 
sition and the final position. For example, to move one 
pseudo- vortex at the plaquette (m, n) in the horizontal 
direction by one site to (m+1, n), we continuously change 
the sign of the coupling parameters J^ +l n , Km+\ n , and 
K^i 4 n . We show in Fig. [4] the spectrum of lowest energy 
levels vs the value of Jf x 10 for moving a pseudo-vortex 
from the plaquette (10, 10) to (10, 11) in an 18 x 18 lat- 
tice. As expected, the spectrum remains almost the same 
after moving the pseudo-vortex. Again, when the size of 
the system and the distance between vortices goes to in- 
finity, the spectrum will remain exactly the same. We 
also give the spectrum versus distance between two vor- 
tices in the same horizontal line by gradually changing 
sequentially the J and K coupling parameters for each 
lattice site, as shown in Fig. [SI The gap between the zero 
fermionic mode and the ground state decreases exponen- 
tially with increasing distance between the two vortices, 
confirming the supposition that decreasing the density 
of vortices leads to an exponential improvement in their 
degeneracy. 

To measure the topological state of the system, we 
propose a method to represent topological information 
by real vortices. Real vortices correspond to negative 
eigenvalues —1 of vortex operator W p = Y[ u jk- I n the 
earlier discussion about the creation and manipulation 
of the vortices, there have only been pseudo-vortices in 
the system: those equivalent vortices realized by chang- 



Figure 4: Spectrum of lowest energy levels vs the value of 
Jfi.io Ior moving a vortex from the top right of the spin at 
(19, 10) to the top right of the spin at (21, 11) in an 18 x 18 
lattice. 




distance 

Figure 5: Excited state spectrum vs distance between two 
vortices in the same horizontal line by gradually changing 
coupling parameters to move one vortex of the pair. Here we 
have set the ground state energy to 0. The horizontal axis is 
the distance between two vortices. The points between integer 
distances correspond to the intermediate states when moving 
vortex from site to site. When the distance is (or 18), there 
are actually no vortices due to periodic boundary conditions. 
Note the exponential suppression of residual energy difference 
as a function of distance, consistent with predictions for the 
topological sector. 



ing certain coupling parameters of Jp?£(t) and K p ^(t). 
Suppose we want to measure the total topological charge 
of two pseudo- vortices, i.e., to measure whether the 
fermionic zero mode shared by the two pseudo- vortices is 
occupied or not. First we move one pseudo- vortex to be 
the nearest neighbor of the other pseudo- vortex. Suppose 
the two pseudo- vortices are at the plaquettes (m — 1 , n) 
and (m, n). At this point, the low-energy fermion modes 
are not degenerate any more-there is a finite splitting 
between the ground state and the first excited state for 
the two neighboring vortices. Next, rather than changing 
J and K together, we first adiabatically change n (t) 
from —1 to 1, with the spectrum shown in Fig. [51 where 



5 



we also give the spectrum for the vortex configuration 
with two pseudo-voritices and two real vortices at the 
same positions. We can see that there the degenerate 
point for the two ground levels of different vortex configu- 
ration occurs at a different value of J than the degenerate 
point for the two first excited levels. 

If, near the degenerate point for the two first excited 
levels, we adiabatically turn on a local term h z 02m- i,n 
by, e.g., turning on a local magnetic field in the z direc- 
tion which acts on the spin at the position (2m — l,n), 
then the two first excited states for the different vortex 
configurations will couple to each other by h z U2m—x,n- 
Denote the first excited state without real vortices by 



\in). As (V>i|of r 



l,n " ™-l,n a 2r, 
2m-l,nl^l) = - 



,_ 1)n |V»l) = -1 and 
1. tffm-i.J^i) has a 
pair of real vortices at the plaquettes (m — 1 , n) and 
(m, n). The local filed results in a finite gap proportional 
to the magnetic field, which we have also calculated nu- 
merically (with a lower bound 0A6h z corresponding to a 
12 x 12 lattice). Note that the gap does not depend on 
the system size if the size is larger than the correlation 
length, since the vortices are next to each other and the 
magnetic field is local. When the system goes through 
the degenerate point for the two first excited levels, the 
first excited state for the configuration without real vor- 
tices becomes the state for the configuration with two 
real vortices. Then we adiabatically turn off the mag- 
netic field. The ground states are not affected as the first 
excited states. 

To finish converting pseudo- vortices to real vortices, we 
finally change K^f n , and K^_ 1 n to complete the fusion 
of the two pseudo- vortices. In this way, after the fusion 
of the two pseudo- vortices, we obtain two real vortices if 
the total topological charge of the pseudo- vortices is 1, 
and obtains no real vortices if the total topological charge 
of the pseudo- vortices is 0. 

We now conclude with a discussion of real- world con- 
straints on topological computing which can finally be an- 
alyzed using this explicit model for braiding and measure- 
ment. First, a honeycomb lattice with periodic boundary 
condition may be difficult to realize experimentally. To 
obtain free boundary condition, we can simply set the 
links connecting the boundaries to 0. 

In order to study the statistical properties of the vor- 
tices, we employ adiabatic theory to calculate the evo- 
lution of the s yste m as we braid the vortices, similar in 
spirit to Ref. [lg. When a system has an almost de- 
generate ground subspace, separated from the rest space 
by a finite gap, and if the Hamiltonian changes adiabati- 
cally compared to the energy gap, the system will evolve 
within the ground subspace throughout the process if it 
is originally in the ground subspace. Let us further make 
the approximation that the ground subspace is degener- 
ate. The state of the system can be denoted by 



*(t) = 2c„(t)Vn(t)expTJ?<0* 



(8) 
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Figure 6: Spectrum for fusion. From t — to 1, the J between 
two adjacents vortices is changed from negative to positive; 
from i = 1 to t = 2, the associated three-spin interactions' Ks 
are tuned from negative to positive. Only adiabatic evolution 
is considered here, corresponding to units of time much longer 
than 1/J. Offsetting these two operations leads to crossing 
between the two sectors (with and without real vortices) at 
different points in time. Addition of a small magnetic field 
during one of the crossings leads to an avoided crossing and 
an efficient means to convert topological information into spin 
information which can be measured through the vortex oper- 
ator W. 



ground subspace. exp"S" /o c (* ) d * is an overall dynamical 
phase factor, while c n (t) contains the geometric informa- 
tion. fy(t) evolves according to 



dc n (t) 

at 



Y.^{t)\ d -^) Cm (t). 



(9) 



By diagonalizing the generalized Kitaev model, one can 
obtain the spectrum and fermion modes of the system. 
The basis states ip n (t) can thus be determined by the 
fermion creation/ annihilation operators. For example, 
when there are four vortices, there are two zero fermion 
modes. The ground state without any fermion excitation 
ipo(t) satisfies 



b k (t)Mt) = 0, 



(10) 



for all k. Vi (t) = b\(t)ip (t), i/j 2 (t) = b\{t)ip Q {t) , and 
if>s(t) — b\(t)b\(t)ipo(t). According to the adiabatic the- 
ory, we can focus on the ground subspace expanded 
by ^0,1,2,3(^)5 and thus &i,2(i) and b\ 2 (t). To calcu- 



late (4>n(t) 



dt 



), we can expand b\^{t + At) and 



where ip n (t) is the nth eigenfunction of H(t) in the 



b\ 2 (t + At) by 61.2(f) and b\ 2 {t). Now we can study 
the statistical properties of the vortices in the general- 
ized Kitaev model. In the following we discuss the main 
results obtained from numerical calculation. 
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The first aim is to prove the SU(2)2 anyonic behavior 
of the vortices. At a first glance of the numerical result, 



e.g. 



B12 

B23 



-0.7155- 0.6969i 
0.0197- 0.0429i 



-0.7141 
-0.3830 



0.05M 
0.5836i 



-0.0197 
-0.7155- 

0.3830 4 
-0.7141- 



- 0.0429i 

- 0.6969i 

0.5836i 
-0.0510i 



(11) 



(12) 



the braiding matrices B12 (by braiding the first two vor- 
tices) and B23 (by braiding the middle two vortices) are 
not like what we expect from SU(2)2 anyon model. Un- 
like the SU(2)2 anyon model, instead of each two anyons 
sharing one fermion mode, the four vortices together 
share two fermionic modes. Even non-braiding relative 
motion between two vortices can induce the evolution of 
the state of the system. In fact, B12 and B23 for the vor- 
tices differ approximately from those for SU(2)2 anyon 
model by a unitary transformation, 



B 12 « C/ f 
B 23 « C/ f 



-e"/ 4 

- e - M / 4 



U 



y/2 sf2 
i 1_ 

\f2 \[2 



u 



(13) 
(14) 



where U = 



0.9576 4 
-0.0335 



0.2859i 0.0335 
-0.0047i 0.9576 



0.0047i 
0.2859i 



(15) 

while the algebraic structure of the braiding operations 
is the same as in the SU(2)2 anyon model. 

For a finite size system, the zero fermionic mode is 
usually not exactly 0, which results in the split of the 
degeneracy of the ground states. Thus there will be dy- 
namical phase errors to the braiding operations. Even 
if the dynamical phase errors are corrected in some way, 
the braiding matrices will still be different from the exact 
form. This finite size effect is usually not discussed, and 
may be used to carry out universal quantum computa- 
tion with higher accuracy. There is another way to make 
the generalized Kitaev model support universal topolog- 
ical quantum computation. We can change the coupling 
parameters suddenly or apply certain sudden magnetic 
pulse on a single spin to transport a vortex. This equiva- 
lently cut the braid in the (2+l)-dimensional spacetime, 
and the evolution differs from the adiabatic evolution. 
The sudden and adiabatic changes of the coupling pa- 
rameters together may also support universal topological 
quantum computation. 

Another interesting result is that there seem to be 
two different topological phases (conjugate to each other) 
for the Kitaev model. In order to identify the different 
phases, we create two vortices and calculate the braid- 
ing matrix B12 (which is a diagonal matrix or a phase 
gate). Choosing Jo = 1, and varying the value of Kq, 
the phase changes from — i to i when K changes across 
some critical value, e.g, between 0.17485 and 0.17490. 
This phase transition is accompanied by an exact zero 




Figure 7: Error induced in the memory due to random varia- 
tion (equally distributed within the ranges %0.2 (blue circles) 
and %0.5 (red stars)) of all coupling parameters. The Hamil- 
tonian returns to the original Hamiltonian after 100 steps of 
random walks. The separation distance between the vortices 
is given on the bottom axes, while the error is the left axis. 
The error is defined to be the distance between the gates G 
and G with and without noise, i.e. d (^G, G^j = \\G — G\\ = 

sup||^|| =1 || (g — G^j ip\\. For the topological memory, G is 
the identity operation, while G is the operation with error. 
Increasing separation distance s between the vortices leads 
to suppression of errors, with apparent exponential behavior, 
consistent with topological protection against local noise. 



fermionic mode at the critical point, where the creation 
operator b\(t) changes to the annihilation operator bi (t), 
and thus ipo(t) and ip±(t) exchange their roles with each 
other. Such critical points form a subset in the parame- 
ter space of the coupling parameters separating different 
phases. 

We can also introduce disorder into the system by 
changing some coupling parameters. For example, we can 
set all the coupling strength of the links around a spin 
to be 0, which creates a vacancy in the system. Our nu- 
merical calculation shows that braiding a vortex around 
a vacancy point is a trivial operation. 

An important question about the Kitaev model is its 
topological stability. It is quite difficult to discuss the 
effect of general noises. We first discuss a class of noises 
which changes the strength of the coupling parameters 
and can be calculated exactly numerically. For example, 
we allow the coupling parameters to deviate randomly 
from the accurate values but within certain range for a 
system with four vortices, which form a single qubit. This 
results in the random walk of the topological quantum 
memory, though there are neither creation/annihilation 
nor braiding of vortices. We can see that when the system 
size is increased, the error of the memory is suppressed 
exponentially, as shown in Fig. [7] We can also allow 
the coupling parameters to deviate randomly from the 
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accurate values but within certain range as we braid the 
vortices, which results in similar robustness behavior of 
the topological quantum gates. 

Combining the above discussion with perturbation the- 
ory, we can discuss the effects of general slow noises which 
act locally on the spins. There are three types of effects. 
One effect preserves the number and positions of vortices, 
but changes the fermion spectrum. Specifically, the two 
and three orders of this effect results in the variation 
of link parameters, which we have numerically studied, 
while higher orders introduce new terms which are prod- 
ucts of the lower order terms. The resulting errors are 
topologically suppressed. The second effect is hopping 
of vortices, which preserves the number of vortices but 
changes their positions. The hopping of vortices can be 
suppressed by introducing six-body interaction (as in the 
case of three-body interaction) which fixes the positions 
of the vortices. The third effect is creation/annihilation 



of vortices. The creation of vortices is suppressed by 
the excitation energy, while the annihilation of vortices 
is suppressed by the distance between vortices. 

To further simplify the generalized Kitaev model, we 
can choose K?,' 3,5 ' 6 (t) = to reduce the number of aux- 

l J j k x ' 

iliary spins. This simplification should be compensated 
by larger system size. 

We have suggested a complete means of analyzing the 
non-Abelian sector of the Kitaev honeycomb model, in- 
cluding the analysis of noise and the observation of topo- 
logical suppression of such noise. Our approach may also 
be applied to other similar systems such as in [19]. How- 
ever, key open questions remain, including understanding 
the effect dynamical (high-frequency) noise plays. This 
would correspond to terms which drive the system above 
the gap, and may be unavoidable in physical systems un- 
der consideration. Nonetheless, such questions may now 
be answered using the framework presented herein. 
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